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GEOS- I IASER RETROREFLECTOR DESIGN AND PRELIMINARY SIGNAL CALCULATIONS 
Introduction 

GEOS will be attitude stabilized along the gradient of the gravitational 
field, so that the normal to the bottom plane surface will point toward the 
center of the earth. On this bottom plane will be mounted an array of cube- 
corner prisms, each with an entrance aperture approximately 1 inch in diame- 
ter. These will form a retro- reflecting surface and will enhance the reflected 
energy received in the vicinity of a ground laser source directed at the 
satellite. The aims of this paper are (a) to discuss the advisability of 
tilting the symmetry axes of the cube- comers at various angles to the plane 
normal in an attempt to equalize the effective reflecting area presented to 
rays incident over a large range of angles, and (b) to perform preliminary 
signal calculations which would be useful for the design of ground stations. 



Variation of Effective Reflecting Area vith Angle of Incidence 


As a satellite moves over the laser reflection station during a pass, 
the angle between the line of sight and the satellite symmetry axis will 
change, as shown in Fig. 1. The individual cube-corner prisms will be 
identical to those used on the S -66 Beacon Explorer Satellite in which the 
corners of the triangular front face have been cut off to form an aperture 
in the form of a regular hexagon. For those prisms, the effective reflecting 
area has previously been calculated by P. 0. Minott (Monthly Research and 
Advanced Technological Development Activity Report for April 1963 ) • In 
addition to the angle of incidence, , the variation of area depends upon 
the width across flats, w, .the depth of the prism, d, (which can be expressed 
in terms of v) and its index of refraction, N. When rotated about an axis 
parallel to one of the sides of the hexagon, the expression for the effective 
area can be written 

A = (w tan 30°) coso6^(w - 2d tanc^) + (w - 2d tano/) 2 sin 30° ^J 

where oL / - arc sin ( — ^ ) (l) 

A plot of the relative projected area versus angle of incidence is presented 
in Fig. 2, in which the area at normal incidence' has been normalized to unity. 
The curve resulting from rotation about a corner of the hexagon instead of a 
side does not vary significantly from this. 

In order to maximize the angle of the cone about zenith in which a 
ground station would find effective laser tracking possible, it has been 
suggested that the GEOS reflector consist of an array of cube corners with 
their axes tilted at a variety of angles with respect to the satellite symmetry 
axis. Thus, if the distribution of tilt angles were chosen properly, it was 
hoped that the satellite might present a reasonably large effective area even 
when the zenith distance, 0, was relatively large. Further consideration of 
the proposal shows that there is little or nothing to be gained by pursuing 
this approach. 


Suppose the zenith distance were such that the line of sight made 
an angle oi with the normal to the bottom surface. Then, if all the cube 
corners were parallel and pointing down, the relative projected area would 
be fj (see Fig. 3)*; Now, if a group of prisms were tilted so as to point 
closer to the line of sight, making an angleO^g with the laser beam, these 
reflectors would have a larger effective relative area fg. At the same 
time, it is clear that another equal group of prisms would have to be tilted 
away from the line of sight by the same angle in order to maintain rotational 
symmetry about the satellite axis . These reflectors would present the smaller 
relative area f^ • Where the function f (oC ) is linear, the loss of area in the 
group tilted away from the laser just offsets the gain experienced by those 
tilted in the favorable direction. Where the second derivative of f(cC ) is 
negative (near zenith), the result is a net loss of reflecting area, and 
where the second derivative of f(v£ ) is positive, (nearer to the horizon) there 
is a net gain. 

The small improvement in effective area that might be experienced at large 
zenith distances by distributing the cube corner directions in this manner 
does not contribute appreciably to the usefulness of the satellite because 
other parameters which also determine the signal strength will be found to 
vary in a much more significant way. In particular, the received' energy falls 
off very rapidly as the satellite moves away from zenith because of increasing 
range and increasing atmospheric attenuation. In addition, any attempt to tilt 
the cube corners without lifting them out of the plane of the bottom surface 
would result in shadowing effects and/or degraded geometrical packing efficiency. 
Our conclusions therefore are to arrange as many reflectors as possible into a 
mosaic on the available area, with all of them pointing parallel to the gravi- 
tational gradient. 


Reflected Signal Vs. Zenith Angle 


We use the following parameters: 


E. = the energy in the transmitted laser pulse 

u 


A c f(* ) 


£ «= half-angle of the transmitted beam divergence cone. We 

assume the energy is uniformly distributed over this cone. 
oC ) = effective retro-reflecting area of the array, for light 
incident at an angle oC to the symmetry axis . A q is the 
maximum projected area and f(o) = 1. Fig. 2 is a plot of 

t(u). 

^ = effective reflectivity 

5 * half-angle of the reflected beam divergence. Again, we 

c 

assume that the energy is uniformly distributed in this 


cone . 

(0) = one-way transmission through the atmosphere for light at 
zenith distance 0. ' 

a = slant range from observer to satellite 
d = receiving telescope aperture 
^ = receiving telescope transmission. 

The ratio of received signal to transmitted pulse energy is then given 


where K = 


is a constant, and 


5. ft' 
pA„ 1 4* 


a KG(e) ; 


0(9) . 


is a variable over the duration of a pass since it is a function of zenith 
distance 0. From Fig. 1 and the law of cosines. 




CU = -R cos 0 + £r 2 cos 2 e + h(2R + h)^ 




From the law of sines, 

, R 

sin oL <= s in 0 


( 5 ) 


( 6 ) 


These equations enable us to calculate expected signal levels for a satellite 
at different positions in the sky. 

Variation of Atmospheric Transmission 

To estimate the effect of zenith angle on atmospheric transmission, we 
assume the very rough model shown in Fig. k. The atmosphere here is a 
uniformly dense layer of thickness T. A light beam passing once through the 
entire atmosphere in a vertical direction will be attenuated so that its 
intensity is reduced to a fraction “2^0 of its initial value. At a zenith 
angle 0, the atmospheric path length is t(0) and the one-way transmission is 






where t(0) a £ -R cos 0 + ^R 2 cos 2 0 + T(2R +T)^j ^ 


( 7 ) 

( 8 ) 


Note that in the limit of a flat earth (R—>°o), equation (6) may be 

written „ 

CV- 2.T" — p-X> 

t(0) “ R cos ejp + ■+* 

= R C-CKL, I •+• ^ ^ \+- — |"£ 


- ' } 


and 


t(0) T sec 0 

.AJLc, & 


( 9 ) 

(io) 


X (e) « 

This is the usual expression used by astronomers for atmospheric transmission.' 
It is only accurate for zenith angles such that 

(u) 


or 


sec 2 © <1 + 0.2 ^ 


If we take for the thickness of our model atmosphere 7" = 25 km, and for 
the radius of the earth R = 6,310 km, we find that the condition in eq. (11 ). 
is satisfied by 

0<8O° (12) 

Since we will never operate at zenith angles even as large as this, we may 
as well use eq. (10) instead of the more complicated expressions (7) and (8). 

The value of is chosen to be 0.7> from several conservative estimates 
in the literature. Pig. 5 is the resulting variation of atmospheric trans- 
mission with zenith angle. 

Example of Numerical Signal Calculations 

The expected orbit for GEOS-I is one with a perigee of 600 nautical miles 
and an apogee of about 1200 nautical miles. We choose an altitude of 1500 km 
(about 800 nautical miles) as a typical example because it results in fairly 
large off-axis incidence angles d, so we can study the drop-off in effective 
area. The corresponding relation between d. and zenith angle 0 is shown in 
Pig. 6. When the effective reflecting area is then plotted as a function of 
0 instead ofet, as required in eq. (4), it appears as in Fig. 7. The last of 
the variable factors, slant-range, Ob, is plotted in Fig. 8. 

We next turn to a calculation of the constant K in eq. (3): The speci- 

fication on the cube corner reflectors is that 50$ of the light incident on 
the effective reflecting area shall be retro-reflected within a cone 20 seconds 
of arc in diameter. We thus have 

f = 0.5 

* 5 

S = 10 sec = 5 x 10 radians. 

The transmitted laser beam will be taken with a divergence equal to that 
planned for S-66 since it seems to be a good compromise between aiming ability 
and energy concentration. 

-4 

» 5 x 10 radians 



The receiving telescope is assumed to have an aperture with a diameter 

d = 16 inches = 0.4064 meters 
and the transmission of the receiving optics is 
^ = 0.5 

The reflector arrays will be fitted into four flat panels which take 
on the trapezoidal shapes shown in Pig. 9- Around the perimeter of each 
panel '/ill probably be a wasted strip about l/2 inch wide because of the 
hexagonal shape of the reflectors . The useable area of these panels adds 
up to about 

A 0 = O.I 78 sq. meters. 

When the appropriate values are inserted into eq. ( 3 ) we obtain 
K > 0.9358 kin*’ # 

The ratio of received to transmitted energy from eq., (2) is plotted in 

Fig. 10. The scale on the right has been converted into received photons 

at the sensor, per joule of. energy in the transmitted ruby laser pulse. 

k 

Using a photomultiplier detector we have found that 10 photons results 
in a good signal,, even against a bright sky, with proper filtering. This 
would allow us, from Pig. 10, to get range information from GEOS down to 
a zenith angle of at least 53° > with much better signals near zenith. For 
other geometrical arrangements, other detectors and lasers, similar calcu- 
lations may be performed. However, it is already clear from the foregoing 
that the laser reflection panels on GEOS will be very valuable for range and 
photographic tracking with geodetic precision. 
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